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Abstract 

We reformulate the one-dimensional complex Ginzburg-Landau equation as a 
fourth order ordinary differential equation in order to find stationary spatially- 
periodic solutions. Using this formalism, we prove the existence and stability of 
stationary modulated-amplitude wave solutions. Approximate analytic expressions 
and a comparison with numerics are given. 
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Introduction 



The cubic complex Ginzburg-Landau equation (CGLe) is a generic amplitude 
equation describing Hopf bifurcation in spatially extended systems, i.e., I Q 
systems [1], with reflection symmetry [6,3,4]. It is of great interest due to its 
genericity and applications to onset of wave pattern-forming instabilities [1] in 
various physical systems such as fluid dynamics, optics, chemistry and biology. 
It exhibits rich dynamics and has become a paradigm for the transition to 
spatio-temporal chaos. 

We consider the one-dimensional CGLe for the complex amplitude field A(x, t): 



A t = fiA + (1 + %a)A xx - (1 + i(3)\A\ 2 A (1) 

where A(x,t) : M 2 i— > C, and fi, a, (3 6 R, x G V. V is the spatial domain on 
which the equation is defined. Interesting domains for us are either the whole 
real axis or a finite box of length L with periodic boundary conditions, /i is the 



Preprint submitted to Elsevier Preprint 



July 29th 2002 



control parameter. Only // > is considered because we study the supercrit- 
ical Ginzburg-Landau equation; one could set fi = 1 by appropriate rescaling 
of the time, space and amplitude, but we keep it as a parameter for closer 
connection with experimental results and previous literature. Coefficients a 
and j3 parametrize the linear and nonlinear dispersion. 

If both a and /3 are set to 0, we recover the real Ginzburg-Landau equation 
(RGLe) in which only the diffusion term and the stabilizing cubic term com- 
pete with each other and the linear term. A Lyapunov functional exists in that 
case [1] and the RGLe behaves like a gradient system. Another limit — the 
nonlinear Schrodinger equation — results from setting cc,/3 — > oo; we then 
have an integrable nonlinear PDE. For parameter values in the intermediate 
range, long-time behavior of the CGLe can vary from stationary to periodic 
and to spatiotemporal chaos [5]. In this paper, we concentrate on the station- 
ary solutions of the CGLe in a finite box of length L with periodic boundary 
conditions, and the case a ^ (3. Stationary solutions are the simplest non- 
trivial solutions, related to propagating solutions by an appropriate change of 
frame of reference (x, t) i— > (x — vt, t) with fixed v e M. 

Searching for coherent structures allows one to reduce a partial differential 
equation into an ordinary one, and such solutions of the CGLe are believed 
to be extremely important in many regimes, including the spatiotemporal 
chaos [9]. Recently, numerical integrations of the CGLe have focused on a 
class of solutions called modulated-amplitude waves (MAWs) and their role 
in the nonlinear evolution of the Eckhaus instability of initially homogeneous 
plane waves [12,13]. 

MAWs can bifurcate from the trivial solution A = (case I) or plane wave 
solutions of zero wavenumber (case II). Analytical aspects of modulated so- 
lutions of the CGLe have been addressed by Newton and Sirovich who have 
applied a perturbation analysis to study the bifurcation in case II [14], and 
discussed the secondary bifurcation of those MAWs [15]. Takac [16] proved 
the existence of MAW solutions using a standard bifurcation analysis in the 
infinite-dimensional phase space of the CGLe, in both cases I and II, together 
with a stability analysis in case I by means of the center manifold theorem. 

In this article we reformulate the CGLe equation assuming a coherent struc- 
ture form for the solutions, and obtain a fourth-order ordinary differential 
equation (ODE) with a consistency condition. This form is algebraically con- 
venient, because the deduced system of four first-order ODEs contains only 
quadratic non-linearity. In the Benjamin- Feir-Newell regime, where plane waves 
solutions are always unstable, we give a proof of existence of MAWs in both 
case I and II using our ODE. For weak perturbations in case I or II, we write 
approximate analytic solutions in the ODE phase space. Coming back to the 
full CGLe, we then prove the stability of those MAWs in a finite box in case 
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II, and prove that the bifurcation is supercritical, as suggested by recent nu- 
merical work [12]. 

In the next section, we discuss symmetries and solutions of the CGLe. In sec- 
tion 2 we transform the steady CGLe for MAWs into an equivalent ODE, and 
give the sufficient condition to identify the solutions of these two equations. 
In section 3 this ODE is used to construct a 4-D dynamical system and prove 
the existence of symmetric stationary solutions of the CGLe in the two cases I 
and II. In section 4 the approximate analytic form of the solutions is given and 
compared to numerical calculations, and the stability of MAWs in case II is 
proved. Several theorems needed in the proofs are reproduced in appendix B. 



1 Basic properties of the CGLe 

1.1 Symmetries 

The equation (1) is invariant under temporal and spatial translations. More- 
over, it is invariant under a global gauge transformation A — > A exp(i</>), where 
(f) e M, and under x — > — x reflection. As a consequence, it preserves parity of 
A, i.e., if A(— x, 0) = ±A(x,0), then A(—x,t) = ±A(x,t) for any later time 
t > 0. If A(x,t) has no parity, then A(—x,t) gives another solution. 

1.2 Stokes solutions and their stability 

The global phase invariance implies that the CGLe has nonlinear plane wave 
solutions of form 

A(x, t) = Ro exp(i(qx — out)) , (2) 

where R$ = /i — q 2 is the amplitude squared, u = fi/3 + (a — (3)q 2 is the 
frequency, and q G K, q 2 < pL is the wavenumber. They are called Stokes 
solutions [2] and are parametrized by the wavenumber q. The two limit cases 
of interest to us are highlighted on figure 1: a plane wave of wavenumber fi 1 ^ 2 
and of vanishing amplitude (case I), and the wave with zero wavenumber and 
maximum amplitude \i (case II). In case II, the solution oscillates uniformly 
in time; we call it the homogeneously oscillating state (HOS). 

For the infinite system, the Benjamin- Feir-Newell [10] criterion states that all 
plane wave solutions are unstable with respect to long wavelength perturba- 
tions (i.e., of wavenumber k — > 0) if 1 + a (3 < 0. If 1 + a (3 > 0, we have to 
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Fig. 1. Marginal stability curve (MS) and Eckhaus instability curve (E) denning 
regions where plane waves solutions exist in the CGLe (inside (MS)), and are stable 
when 1 + a(3 > (inside (E)). Case I corresponds to the wave of maximal possible 
wavenumber and Case II to the wave with q = and the maximal amplitude. 

consider the Eckhaus instability criterion; only a band of wavenumbers are 
stable against long wavelength perturbations (figure 1): 



?<& = ^%- (3) 
3 + ap + 2p 2 



For a finite periodic system the wavenumbers for both the original states 
and the perturbations are quantized. These criteria have been reexamined by 
Matkowsky and Volpert using linear stability analysis [18]. 



1.3 Coherent structures and MAWs 



Coherent structures play a very important role in the study of pattern forma- 
tion and dynamical properties of the CGLe [9]. They are uniformly propagat- 
ing structures of the form 

A(x, t) = R(x - vt) exp(i(f)(x - vt)e~ iuJt ) 

which can be expressed as solutions of a 3-D nonlinear dynamical system 
obtained by substituting the above ansatz into the CGLe. There are two free 
parameters: the frequency to and the group velocity v. 

The fixed points of the 3-D system are the plane waves described in the 
previous section. The homoclinic [11] and heteroclinic [9] connections between 
the fixed points correspond to localized coherent structures. The Nozaki-Bekki 
solutions [7] belong to this category; they connect asymptotic plane waves 
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with different wavenumbers. In numerical simulations in large domains, nearly 
coherent structures are frequently observed in chaotic regimes, thus suggesting 
those objects are also relevant to spatiotemporally chaotic dynamics. 

Recent numerical studies reveal another kind of coherent structure: modulated 
amplitude waves (MAWs) for the CGLe [12]. They correspond to limit cycles 
of the 3-D nonlinear system. When v — 0, MAWs are stationary. The for- 
mation of MAWs is the first instability encountered when a plane wave state 
crosses the Eckhaus or Benjamin- Feir stability line. The MAW structure is 
frequently observed in experiments [6,20] and considered as a key to interpre- 
tation of patterns and bifurcations exhibited during the system's transition to 
spatio-temporal chaos [13]. Traveling MAWs have been observed in numerical 
simulations of the CGLe in periodic boxes, with parameter q between and 
/i 1 / 2 , i.e., in between cases I and II; we are interested here only in stationary 
MAWs that appear either in case I or case II. 

In this paper, we propose a new real-valued ODE to describe steady solutions 
of the CGLe. A A-D dynamical system derived from this ODE enables us 
to apply the successive approximation method [8], to prove the existence of 
stationary MAWs and to give the analytical form of the approximate solutions 
in both case I and case II. Numerical integrations of the exact CGLe are 
then compared to the approximate analytic result. Furthermore, we show non- 
analyticity at discrete points of solutions in case I, and prove the stability of 
the MAWs in case II. Some theorems needed in our proof are reproduced in 
the appendix B. In what follows, diag(- • • ) denotes a (block) diagonal matrix 
and col(- • • ) a column vector. 



2 Stationary case 

Since we are only interested in the steady solutions of the CGLe, we substitute 
the ansatz 

A(x, t) = R(x) exp(i(f)(x) - iut), (R, 0) G K 2 (4) 
into (1). We then have 



(1 + a 2 )G x = K=(/3- a)R 4 - (u - fia)R 2 (5) 
(l + a 2 )G 2 =M = (1 + a 2 )R 3 R xx + (aw + fi)R 4 - (1 + a(3)R 6 . (6) 

where G = <p x R 2 is reminiscent of "angular momentum". Note that if a — /3, 
this "angular momentum" is conserved — it is constant in space — provided 
that us = jia. In that case, (6) can be solved in terms of elliptic functions [17]. 
We will only consider the case a ^ (3 in the following. Equations (5) and (6) 
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are invariant under (G, x) — > (— G, —x). Note that for plane waves, K — and 
G is a constant. If if is not always zero, differentiating (6) and dividing the 
result by (5) gives 

2G = M X /K, (7) 

and by (6) 



l + a 2 M 2 

M =—lt-- (8) 



Furthermore, we can factorize R 2 from M x and K and write M x = R 2 N and 
K = R 2 P, where 

N = (l + a 2 ) l -(R 2 ) xxx + (auj + /U )2( J R 2 ) X - (1 + a(3)?>R 2 {R 2 ) x 

P = ((3-a)R 2 -(to- fwi). (9) 

The last relation can be used to express R 2 in terms of P: 



g= ( "7 a)+p =« + rfp=flg + /-, do) 

p — a p — a 

where d = — a) and e = (cu — /ia) / (0 — a) . 

Note that e = i?o is the square of the homogeneous amplitude -Ro(<?, u;) of the 
Stokes plane wave solution (2) of frequency u and wavevector q(u). P then 
appears as the modulation of the amplitude squared with respect to the Stokes 
solution, and so it is an appropriate variable to describe a MAW. 

Substituting K and M x into (8), we have 

1 + a 2 N 2 



4 P 2 



= M. (11) 



If P 7^ (11) is equivalent to (5) and (6). It is easy to check that if we regard 
(7) as a definition of G, and use K, M, N, P expressed in terms of R, equation 
(5) and (6) will be recovered as a result of (8) and (11). Differentiating both 
sides of (11) results in 



a 2 



(PN X -NP X ) = R 2 P 3 . (12) 



In this step we have extended the solution set of (11), because as we integrate 
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(12) back, we get 



1 + a 2 N 2 

where C is an integration constant. Only when C = 0, a solution of (12) is a 
solution of (11). For this reason, when obtaining solutions of (12), we have to 
check the consistency condition 

^^-^ = (14) 

to make sure that we have a solution of (11), thus a solution of (5) and (6). 
Note that if K vanishes we have to go back to (5) and (6), since in that case 
(11) is not well defined. Let us rewrite iV in terms of P: 



N = {aP xxx + bP x + cPP x ) , (15) 

1 + cr 

where a, b, c are constants 



_(l + a 2 ) 2 
a ~ A((3-a) 

1 + a 2 f2(ocuj + fi) 
2 V 

3(1 +a(3)(l + a 2 ) 
C= 2(/3-a) 2 

After some algebra (here relegated to appendix A), we get an equation for P 
only: 




= -M + kP, M = XP XX + dP 2 + eP . (17) 



3(1 + a(5){uj — jj,a) 



(16) 



A is a fixed real constant that depends on a and (3 only, and that takes two 
different values given in appendix A. A is a transient variable used in the proof 
and derivation but our solutions to the CGLe do not depend on A and do not 
distinguish the two values of A (see section 4). e and k are real parameters 
introduced as e + = e. So (17) has two free parameters: u, introduced by 
the ansatz (4) as the carrier frequency of the solution, and k. the consistency 
condition (14) fixes one parameter. 
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3 4-D dynamical system and the existence of periodic solutions 



Let us take r as the spatial variable, P = P(t) in (9), and rewrite (17) as a 
system of first order equations in r. With N = M T /P and Q = P T , from (17) 
we have 



' M = 


NP 


N = 


±M + kP 

a 


P = 


Q 




\{M-dP 2 -eP) 



(18) 



where the dot represents the derivation with respect to the spatial variable r. 

It is easy to check that P = is a solution of the original equations (5) and 
(6), corresponding to the plane wave solution of the CGLe with frequency uo. 
We will study the behavior near P = and prove the existence of periodic 
solutions for small P. In the CGLe, this corresponds to a weakly modulated 
amplitude wave which bifurcates from a plane wave solution. If P ~ e, where 
e is a small parameter, so are M, N, Q by their definitions. Write 

(M, N, P, Q) = (ex, ey, ez, ew) 

and set k — ki + ek 2 . Substituting these into the 4-D system, we have 

/ 



d_ 



(x\ 




fx) 




1 yz ^ 


y 


= A 


y 




k 2 z 


+ e 




z 




z 







[wj 




[wj 







where A 



0^ 

a 1 

1 



io^o 



The linear part A describes the behavior of the system in the neighborhood 
of the trivial fixed point (0,0,0,0). Note that the system is invariant under 
(t,y,w) — > —(t,y,w). We use this property to simplify our analysis. More- 

— * 

over, this system defines an incompressible flow since V • X = 0, where 
X = (x,y, z,w). It follows from (13) that the system has one integration 
constant C. This constant induces a foliation of the phase space into three- 
dimensional manifolds. Physical solutions, i.e., the solutions of the original 
CGLe, are restricted to C — 0, the manifold that satisfies the consistency 
condition (14). 

These properties strongly restrict the possible distribution of eigenvalues of 
A. We restrict our analysis to the case e/A > 0, then A has eigenvalues 
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{0, 0, icui, — iui} with u\ = \Je/\. In that case, periodic solutions or MAWs 
can exist as we will prove in the following. The evolution of the system along 
either of the two degenerate eigenvalue directions respects the constant C 
foliation: if the solution is on a constant C manifold at initial time, it remains 
there for any later time. 



We now discuss the condition e/A > in terms of an instability of the un- 
derlying plane wave. We can rewrite e/A using (16) and (A. 2). Assuming that 
the solution we are searching for is close to a plane wave, we can use the 
wavenumber q instead of the frequency u, using the dispersion relation (2) for 
plane waves: 



e b 



A a 1 + a 2 
2 



1 + a 2 



If we write 



we have 



2(au + ji) 



3(1 + af3){u) — ixa) 
f3 — a 



3 + a[3 + 2a 2 )q 2 - (1 + a/3)fi 



<1m 



(1 + af3)n 
3 + a(3 + 2a 2 



(19) 



(20) 



- > 0^ 
A 



Vg e [-01, ^fji] 



q 2 >q^ if (l + a/3)>0 

if (l + a(3) < -2(1 + a 2 ) 
if -2(1 + a 2 ) < 1 +a(3 < 



<f < Qm 



(21) 



The corresponding regions are illustrated on Fig. 2. Note that qM.(o(, (3, /i) = 
q^iP, O-i If = \/3\, the positivity of e/A is assured when the corresponding 
plane wave is Eckhaus unstable. If |a| ^ the positivity does not coincide 
anymore with the Eckhaus criterion; this is not surprising considering that we 
do not restrict our analysis to long wavelength perturbations of plane waves, 
but that the solutions we are seeking may have any wavenumber. 

In the following we distinguish two cases. In the first case eigenvalue has 
a simple elementary divisor, i.e., has two distinct eigenvectors citejhale. This 
coincides with case I: the MAW solution bifurcates from the A = state, 
with u ~ fia and hence e/A ~ Afi > 0, for fi > 0. In the second case, 
eigenvalue has only one eigenvector. This coincides with case II: the MAW 
is superimposed over a plane wave with uj ~ so q ~ 0, and 



A 1 + a 



9 



Fig. 2. Left: wavenumber distribution of stationary MAWs in the (a, (3) plane. In 
(BFS), MAWs exist if q 2 > q^. In (Ml), MAWs exist Vg. In (M2), MAWs ex- 
ist if q 2 < q^. Right: regions of existence of MAWs in the (q, fj) plane in the 
Benjamin- Feir-Newell stable regime ((BFS) region). (MS) is the marginal stabil- 
ity curve, (E) is the Eckhaus instability curve and (M) is existence curve defined 
by (20). Stationary MAWs exist outside (M). 



The positivity is insured if the system is Benjamin- Feir-Newell unstable, (1 + 
a/3) < 0. 

In terms of M, N, P, Q, the consistency condition (14) can be written as 



a ~ ^ 2 



1 + cr)M = ( -N - \Q ) (22) 



where in new variables 



M= d{l + 0?) {dp + e)(j - _ i p,_.p ) _ ^(1 + O 2 ) q2 

(acu + fx)(dP + e) 2 - (1 + a(3)(dP + e) 3 . 
Recalling (6), we may express G by 

1 + or 

Here we are allowed to fix the sign of the right hand side expression because 
of the (G,x) i— > (-G, —x) reflection symmetry of (5) and (6). 
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3.1 Case I 



We want the eigenvalue to have non-degenerate eigenvectors, for this, we set 



Consequently, we have 



e = e + -fc = — k 2 . (24) 
A A 



Notice that e ~ to the zeroth order, so R ~ and uj ~ //a, which means that 
the solution to be considered bifurcates from the zero solution A = 0, corre- 
sponding to a plane wave around the marginal stability curve, with wavenum- 
ber q ~ i/i 1 / 2 . This solution is therefore outside the Eckhaus stability region 
when 1 + af3 > 0. 



Let 



(~e\ 



The four eigenvectors of A are: 

(*\ 
1 






1 



( o N 

— ikiUJi 1 



V 



1 



/ 



/ 







D = 



e 
100^ 

a 

110 
1 



D~ 



1 



ikitu 1 1 
1 

a- l \ 2 



V 



e- 1 



- 1 1 
1 



/ 



and X = (x, y, z,w) = D 1 X. The dynamical equations for the new variables 
become 



d7 



-i(x + £w)(y + z) 
-fiy + zf 
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where 



M (wi) = D- l AD = 



^ 0^ 

1 

^ -u\ J 



The angular frequency of the solution Q should be close to u>i, Q 2 = uj\ + ej, 
with the shift 7 to be determined later. Next, we change variables to: 



X — X\ 

V = x 2 

z — Zi sin Qt + z 2 cos Qr 
w = Qzi cos Qr — Qz 2 sin Qr 



(25) 



The 4-D system of equations then takes form: 



X\ = e 



X2 



Z2 



k 2 (x 2 + zi sin Or + z 2 cos Or ) + — (x 2 + zi sin Or + z 2 cos Or ) z 

a 



e 
e 
e 
O 



nX \ r. 

xi H (zi cos Sir — Z2 sinSZr) 

a 



(X2 + Zi sin Or + Z2 cos Or) 



(x 2 + zi sin Or + z 2 cos Or) 2 cos Or 



+j(zi sin fir + z 2 cos fir) cos Or 



o 



oa 



— — ( xi H (^i cos Or — Z2 sin fir) (X2 + 21 sin Or + 22 cos Or) sin Sir 



e 

o 



(X2 + Zi sin Or + z 2 cos Or) 2 sin fir 



(26) 



—7(21 sin fir + Z2 cos Or) sin Or 







OA 



— — I xi H (^i cos Or — z 2 sin Or) (X2 + Z\ sin Or + z 2 cos Or) cos Or 



The proof of the existence of weak MAWs close to P = relies on a series of 
theorems from J. Hale's monograph [8]. We reproduce the relevant theorems 
in appendix B, and refer to them as the need arises. 

Note that the transformation (r, xi, x 2 , zi, z 2 ) — > (— r, — xi, x 2 , — zi, z 2 ) leaves 
the system (26) invariant. So, by definition B.l of appendix B the system has 
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the property E with respect to Q, with 

Q = diag(-1, 1,-1,1). 



As we are interested only in the solutions with definite parity, we may start 
the iteration with the vector 

— # 

X = (0, a 2 , 0, a 4 ) . 

According to Theorem B.4, our solution z(t,X ,€) has the property 

Qz(-t, X , e) = z(t, X q , e) , 

which means that our solutions are either symmetric or antisymmetric. Ac- 
cording to Theorem B.5, the second and the fourth determining equations are 
always zero for this starting vector. For the first and the third determining 

equations, the zeroth order solution of X , i.e. X , may be substituted, and 
we get 



k 2 a 2 + — ( a 2 2 + ^al) = (27) 

CI Zi 

7 \ 2 tta 2 a 4 d 

2n ai + ^a^-\n a2a * = - (28) 

From (28), we have two possibilities: either 04 = or 

f\ 2 n 2 2d\ n 

7 + «2 — --^ =0. (29) 
V ae a I 



When 04 = 0, using X = (0, a 2 , 0, 0) in (26) leads to a trivial constant solu- 
tion. In the following, we consider only the second case (29). We can solve (27) 
and (29) for 7 and a 4 and prove that the system (26) has periodic solutions. 
Note that we have three free parameters e, a 2 ,k 2 . But as we will see further, e 
and a 2 are always combined as ea 2 in the first approximation controlling the 
amplitude and the period of the solution, and the combination will therefore 
be regarded here as one single free parameter. For general periodic solutions, 
a 2 can be interpreted as a phase control parameter, i.e., a parameter giving 
the initial location on the periodic orbit at r = r . Here, because we only 
consider symmetric solutions, the translational symmetry of the autonomous 
system is broken, and that is the reason why e and a 2 combine into a single 
parameter. The remaining parameter k 2 can be chosen freely, for example as to 
satisfy the consistency condition (22), which, when the zeroth order solution 
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is substituted, becomes at order (e 2 ): 



-*Ltfal + ii(da 2 + =0. (30) 



At zeroth order, Q 2 = uj\ = 4/x and e = A/i\. Solving the system of equations 
(27), (29) and (30), we get 



k 2 = - — da 2 

7 : 

CLa 



a°2 

±2a 2 



We can write out the Jacobian for those three equations explicitly: 

/ 

J = 



a 2 






1 



— a a 

a * 



2f (da 2 + *f) -f ^ 2 a 4 



fc2a 







The determinant of this Jacobian is 

1 



det J = -d 2 Q 2 a 2 aA ^ a 2 ^ . 
2 



We now invoke theorem B.2, reproduced in appendix B, and conclude our 
proof that system (5) and (6) has periodic solutions near P = 0. We shall give 
approximate solutions in section 4, and show that in this case they contain 
defects. 



3.2 Case II 



Eigenvalue has only one eigenvector. In this case, we assume that -e + k x ^ 
to the zeroth order in e, so without loss of generality we can choose k 2 = 0. 

Then = + k\. Implementing the transformation X = DX with 



e 





1 


—kiUJi 


1 1 








1 
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we have 



where 



^-X = M(coi)X + e 
dr 



\(x - kihw)(y + z) 
-\{x - h±w)(y + z) 



M(wi) = D~ l AD 



/ \ 

a 


1 

o o -u\ o y 



As in case I, let Q 2 



cut 



C7 and perform the same transformation (25) into 



variables Xi, x 2 , Zi, z 2 . We then obtain a A-D system similar to (26). However, 
in the equation for x\, there is an e-free term. In order to use the succes- 
sive approximation method, further transformations are required. Let p G H. 
such that p 2 = e. With the transformation x 2 — > px 2 ,e — > p 2 we recover the 
standard form 



Xl 
X 2 

Zl 



a 

P 



x 2 — p 2 k\ — (px2 + Zi sin fir + z 2 cos Or) 2 



oa 



xi — k\—z-{z\ cos Or — z 2 sin Or) 



(px 2 + zi sin Or + 22 cos Or) 



El 
o 



-d 

T 



(/)X2 + zi sin Or + z 2 cos Or) 2 cos Or 



+7(21 sin Or + z 2 cos Or) cos Or 



(31) 



(32) 



— - (xi —ki (zi cos Or — z 2 sin Or) i (/; v 2 + : 1 sin Or 4 z 2 cos Or) sin Or 



z 2 = 7T 







v (p^2 + ^1 sin Or + z 2 cos Or) 2 sin Or 
A 



—7(21 sin Or + z 2 cos Or) sin Or 

— — (yXi —ki [z\ cos Or — z 2 sin Or )^ (px 2 + zi sin Or + z 2 cos Or ) cos Or , 

(33) 

The system (32) has the same symmetry as identified in the case I. If we 
are only interested in solutions with definite parity, we may again start the 
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iteration with X = (0,a 2 ,0,a 4 ). To the second order (p 2 ), the determining 
equations are: 



a 2 ^- P ^ + 0(p 3 )=0 (34) 
7a 4 2 (da 2 a A A^a 2 a 4 /ci\ 3 

From the second equation we obtain either a 4 = (trivial for our purposes, 
as discussed above) or 

(2d XWkA , 
7 - P«2 f y + -^r 2 ) + 0(p 2 ) = . (36) 

If we backtrack the transformations made, it is clear that the consistency 
condition requires that we keep terms up to the fourth order (p 4 ). We found 
that with the substitution 

auo + p ,,2 . 
e = — — - + p (p w 3 - pda 2 ) , 
1 + ap 

where is a new parameter, only the fourth or higher order terms are left 
in the consistency condition. From the definition e = Rq = {u — /ia)/((3 — a) 
and the above equation, we get u ~ p/3 and then e ~ p to the zeroth order. 
So i?o ~ y^) Q ~ 0) which means that this solution bifurcates from the 
HOS A = v /pexp(— iujt). To the leading order (p 4 ), we are allowed to use the 
following substitutions in the consistency condition (22): 



a 2 — > c<j — > p/9 f2 



2p(l + 



1 + a 



2 



gW 2A(l + «fl)\ 2pA(l + a g) 

a 1 + a 2 / 1 + a 2 



The resulting equation is of a relatively simple form: 

a\(-X + d 2 (l + a/3)(l + a 2 + \ + \a(3)) + 4(1 + a/3) 2 Apw 3 = . (38) 

From (34) it follows that a 2 is of order p, and from (36) that 7 ~ 0(p 2 ). After 
a change of variable a 2 = pa 22 and keeping only the highest order for the 
equations, we can rewrite (34) and (36) as 



Ae hdaj 

a 2 2 ^^ = ( 39 ) 

a 2e 

7 = 0. (40) 
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For e,e,ki we use the values in (37). From (39), (40) and (38), we can solve 
for a 2 2,7,^3- The Jacobian of those equations is 








\ 



J 



1 







V 



4(1 + a(3) 2 \fi 



) 



So, det J = 4(1 + a(3) 2 X 2 fie/a ^ for 1 + a(3 ^ 0. According to Theorem B.2, 
we have proved that equations (5) and (6) possess periodic solutions. 



4 Analytic form of periodic solutions, stability analysis and nu- 
merical tests 

We have proved in the preceding section the existence of symmetric periodic 
solutions in case I and II. In both small parameter e or p ensures the 

convergence of successive approximations. However, we did not give a bound 
on the highest value of this parameter, nor did we show that the solutions 
which we obtain are the ones observed in numerical simulations. In this section 
we give the approximate analytical form of periodic solutions. We compare 
them with direct numerical integration of the CGLe in case II. 

The solutions are shown to be independent of A to order e in case I and to 
order e 2 in case II. In addition, these solutions should also satisfy the 3-D ODE 
mentioned in section 1.3 which do not contain A, so they can be matched with 
the solutions of the 3-D system in a unique way, independent of the value of 
A. Hence, we conclude that to all orders the physical solutions are identical 
for the two values of A. 

The two cases are taken separately. In this section, we reinstate x as the spatial 
variable, R = R(x). 

4-1 Case I 

Using (10), (23) and the case I calculations of the preceding section, we have 
after some algebra: 




(41) 



<j>x = 



ea,2 sin 2Vtx ± 2 sin Vtx 



2(l + a 2 )tt 1± cos VLx 
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To the first order of e, R and 4> x are independent of A. The ± sign selects 
two solutions which transform into each other by translating by a half period. 
This is reminiscent of the spatial translational invariance in the symmetric 
solution space. From the definitions of e, f2 and from (24), (19), we get to the 
first order: 



LU = jJLOi 



4/x + 



Q = io\ + 



3ea2 
Geda,2 
1 + a 2 

2U)! ' 



ap + 2a 2 + 3) 



(42) 



We see that u and Q are independent of A. On the other hand, for periodic 
boundary conditions, we can use Fourier modes directly to transform the PDE 
(1) to a finite set of approximate ODE's by Galerkin truncation. Then the 
stationary solution can be obtained by solving a set of nonlinear algebraic 
equations. 



Numerical comparison 

If we take as an example the following parameter values (previously used 
in [21]) for which defect chaos is expected: 

a = 1.5, (3 = -1.2 

and fix the size of the domain to L = 24, then at fi = 0.072644, oj = 0.097879, 
a periodic solution of period L/2 is found. This solution has R max — 0.0750. 
On the other hand, if we use the same a, (3, \x and search for R max ~ 0.075 by 
adjusting e (we always keep a 2 = 1), we find that 

e ~ 0.00380, u = 0.097566, period ^ = 12.0102 . 



The approximate analytic solution and the numerical solution of the exact 
CGLe agree very well. The profile of R from our successive approximation is 
shown in Fig. 3. 



Structure near the defect 

It is easy to see from (41) that only ea 2 > is the physically interesting 
combination. However, we may wonder whether it is really true that R 2 = 
dP + e remains non-negative everywhere while touching zero at some points. 
Fig. 3 and the first equation of (41) suggest a positive answer to this question. 
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R 




Fig. 3. Spatial profile of the amplitude R(x) at u = 0.072644, from (41) with 
Rmax = 0.075. 

But since we have only an approximate solution, further justification is needed. 
Suppose at some instant x , we have dP + e = on the periodic orbit. From 
the consistency condition (22), at this transition point 



d 2 (l + a 2 ) 



a ~ 



Q2+^1 N -XQ) = 0, 



so, Q = N = 0. According to (18), M = and P = 0. Assume that Q = 0, 

then N ^ since the point is not an equilibrium. At next instant xq + 8x, the 
consistency condition can not be satisfied as the two sides of (22) have different 
orders of 8x. So we conclude that Q ^ at the point Xq, which means that 
Q(xq + 5x) has negative sign to that of Q(xq — 5x). Thus, after touching the 
zero value plane, dP + e returns to the positive half space again. The turning 
happens exactly on the dP + e = plane. We claim that dP + e > always 
holds and the equality holds periodically. From (41), in the neighborhood of 
R = at x = xq on the periodic orbit, R behaves like 



R~ 




1/2 



\X — Xq\ 



and is manifestly a non-analytic function of x. 



We do not discuss the stability of the solutions in case I, as this has already 
been accomplished by Takac [16] who has proven that these solutions are 
unstable. 
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4.2 Case II 



To the first order of e, the solutions are 

Xl = - e ±^{2de + AA; 3 )sin2fi:r 
x 2 = t{a\ — A ^° 4 cos 2 fix) 

z x = - 12g £ 2 a ^ 2 (3(3di 2 + A 2 fi 2 fci) sin fix + (rfg 2 - A 2 fi 2 fci) sin 3fix) 
z 2 = a 4 + i 2 g°AH^ (3(A 2 fi 2 /ci — <ie 2 ) cos fix + (A 2 fi 2 /ci — <ie 2 ) cos3fix) , 

where e = p 2 > 0, and 04 is a free parameter. In the following, we will see that 
e and 04 always emerge in the combination e a 4 . To the second order, u is 

e 2 a\ ( 1 + a 8 3 - a \ 
u = ^ + 1 n ; — 9T -5 + 



4p(l + a 2 )\ v /3-a 1 + adJ ' 

It is independent of A, and therefore e, 6, fi are also independent of A. R and 
4> x can also be calculated to the second order: 



d 2 

R 2 = (ec^) 2 + rfea 4 cos fix + 



2/i 



12fi 2 



- + % j cos2fix + e (43) 
a 0/ 



02 



ea4 



pfi(l + a 2 ) 



e sin fix 



ea 4 / „ , o2 7e 



24fi 2 



2 7ce\ 

6dfi^ + -— + sin 2fix 

afi" 1 a J 



So clearly i? and (j) x are independent of A. Similarly, the different signs of <x\ 
will give the same solution up to a half-period translation. This solution is 
the one observed in the numerics when passing the Eckhaus instability for 
underlying wavevector q — 0. Linear stability analysis reveals [18] that the 
q = state, the most stable state under the long wavelength perturbations, 
becomes unstable when the size of the system is such that the smallest possible 
nonzero wavenumber k satisfies 

2 2p(l + af3) _ 2 

K < — - = K . 

1 + cr 

It is easy to see that k 2 = uf up to order (p 4 ). 

For our parameter choices fx = 1, a = 1.5, = —1.2, the bifurcation size of the 
system is L Q = — = 8.95492. In the following, we will first prove the stability 
of our solutions near the bifurcation point. Then we will compare them with 
the stable solutions observed in numerics. 
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Stability analysis: presentation 



Assume that A = Rexp(i(p) where R, G R is an exact solution of (1). The 
perturbed solution is assumed to be A = (R + r) exp(0 + 9), where r, 9 G K 
is the perturbation on the amplitude and phase, separately. Substitute it into 
(1), keeping only the linear terms in r and 9. We have 



r t = (n - - a<p xx - 3R 2 )r + r xx - 2a<p x r x 

-(2R<j) x + 2aR x )9 x - aR9 xx (44) 
R9 t = (u - a4> 2 x + XX - 3f3R 2 )r + ar„ + 20 x r x 

+ (2/4 - 2aR4> x )9 x + i?0 xx , (45) 

where in (45) we have used <p t = —uo. To study the stability of the starting 
solution A, we treat these equations as an eigenvalue problem for a two com- 
ponents vector, i.e., we let r t = o~r, 9 t = o9 and we investigate the spectra a 
of the linear operator resulting from (44) and (45) in the C 1 continuous peri- 
odic function space. As the CGLe has global phase invariance, the eigenvalue 
equations always have solution (r, 9) = (0, 9 ) with eigenvalue a = 0. At the 
same time, spatial translational invariance implies that another eigenmode has 
o = 0. As a result, saying that the solution is stable means that it is stable 
up to a phase and a spatial translation, and that all other eigenmodes have 
eigenvalues with negative real parts. 

Invoking the expression for R, <p x to the second order of e, the coefficients of 
various terms of r, 9 and their derivatives in (44) and (45) become explicit 
functions of x. The resulting linear operator on (r, 9) has even parity due to 
the symmetry of our solution, and we can consider the even and odd solutions 
of r,9 separately. If we set e = 0, i.e., the starting state A is a plane wave 
state, then cos(nQx) and sin(nflx) are the eigenfunctions of the unperturbed 
linear operator. They give the stability spectrum of the plane waves. Now, let 
us move a little (to the order of e) beyond the bifurcation point. The eigen- 
functions are still cos(nflx) and sm(nflx) up to e corrections. For example, if 
the even solutions are considered first, we assume that to the first order the 
eigenfunctions are (the time dependence for r, 9 has been suppressed): 



r = m\ cos(nf2x) + e(m cos((n — l)Qx) + m 2 cos((n + l)Qx)) (46) 
9 = n\ cos(nfirr) + e(n cos((n — l)Qx)) + n 2 cos((n + l)Qx)) , (47) 

where n is a non-negative integer. Note that we do not include the terms such 
as e 2 cos((n±2)Qx) in the above expressions because they induce corrections of 
order e 3 or higher in the eigenvalues. Now if we substitute (46) and (47) into 
the eigenvalue equations and identify the coefficients of cos(nf2x), cos((n — 
l)Qx) and cos((n + l)Qx), a set of six homogeneous linear equations for 
m , m 1 , m 2 , n , ni, n 2 can be derived. The determinant of the coefficient ma- 
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trix will give an eigenvalue equation for a. The resulting expression is too 
complicated to merit being displayed here. 

Before bifurcation, the HOS is stable. The first instability occurs for n = 1 
mode, one eigenvalue of which is very close to near the bifurcation point, 
being negative before and positive after. Meanwhile, for n > 1 modes, the 
corresponding eigenvalues have negative real parts bounded away from zero. 
As the bifurcating solution emerges continuously from the HOS, near the 
bifurcation point (e<l) the perturbed linear operator has all the eigenvalues 
with negative real parts away from for n > 1 and one eigenvalue close to 
for n = 1. So, we only need to check the stability of our solutions for n — 1. 

For convenience, we can fix parameters a and (3 to any values allowed by (21) 
and perform the above stability analysis of the solution. 

Stability analysis: numerical checks 

The numerical values we used are /i = 1.0, a = 1.5, (5 = — 1.2,a 4 = 1. The 
eigenvalue equation is then 



(j = corresponds to the neutral mode associated with the global phase 
invariance. All others solutions have negative real parts. The a_ = — 0.14154e 2 
solution is the interesting one. If we use the same parameter values to calculate 
the stability of the HOS, the eigenvalue equation for n = 1 is 



To the second order in e, we have a = -2.98462 -0.19325e 2 or a+ = 0.07077e 2 . 
The later positive eigenvalue indicates that the plane wave solution is not sta- 
ble. We note that 2a + = —a_ to order e 2 which indicates a supercritical 
pitchfork bifurcation. We have proved that this equality holds exactly at the 
bifurcation point for any values of a and j3, and this justifies the above nu- 
merical checks. Under perturbation the HOS will evolve to the modulated 
amplitude solution given above. When the instability is saturated, the corre- 
sponding eigenvalue for the MAW is negative. If we change the sign of a 4 or 
use the other value of A, the eigenvalue does not change, as expected. 

If we alternatively consider the odd-parity function space {sin(nf2a;)} ne N, we 
obtain the following eigenvalue equation: 



7.9860eV + (56.423 - 63.394e 2 y 2 + (82.564 
+ (45.022 - 28.859e 2 )a 4 + (10.923 

+ (1.0- 



- 75.135eV 3 

- 4.3059eV 5 
0.18864e 2 )a 6 = 0. 



e 2 (-0.21122 - 0.26402(x) + 2.98462a + a 2 = . 
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(28.2115 - 33.6568e 2 V + (27.1763 - 21.2703e 2 )<7 2 + (8.92308 - 

4.04125e 2 )cr 3 + (1 - 0.19287e 2 ) ( r 4 = . 

This equation is quartic because for n — 1 only two modes sin fix and sin 2Qx 
are used. Now a = corresponds to the neutral mode associated with the spa- 
tial translation of the CGLe. Other eigenvalues of the equation have negative 
real parts bounded away from 0. 

To summarize, our solution is stable in the whole phase space of the CGLe, 
up to a phase and a spatial translation. 

In ref [6], B. Janiaud et al. have investigated the stability of traveling waves 
near the Eckhaus instability in Benjamin- Feir stable regime. They derived 
a necessary condition for the bifurcation to be supercritical and located the 
corresponding regions as two strips in the a, ft parameter space. We have 
studied the stationary MAWs in the Benjamin- Feir unstable regime and found 
that the bifurcation from the HOS to MAWs is always supercritical, even when 
parameter values lay outside of the region given in ref [6]. 

In ref [15], application of the perturbation method to the zeroth order (e°) 
equation gave nonzero eigenvalue Xq = 2//3. This can not be correct since the 
zeroth order equation just gives the stability of the unstable HOS. Further- 
more, in the Galerkin projection, somewhat surprisingly the N — 1 truncation 
was found to give a better result than the N = 2 truncation. In our case, if we 
use only the first order expressions for R,(f> x in (44) and (45), we cannot get 
the correct eigenvalues even near the bifurcation point, not to mention that it 
would not be possible to extend the result to the next bifurcation. 

Comparaison with numerical integration of the CGLe 

In our numerical simulations we employed a pseudo-spectral method to evolve 
equation (1) using 128 modes. For system size L < L , we always recover the 
HOS (q — 0). For L slightly larger than L , however, the solution relaxes to 
the modulated amplitude solution given irrespective of the initial condition. 
Figure 4 depicts the stable steady solutions given by the two methods. 



5 Conclusion 

We have reformulated the stationary one-dimensional CGLe in a finite box 
with periodic boundary conditions as a fourth-order ODE for a variable P that 
can be interpreted as the modulation of the amplitude squared of a plane wave 
solution. This reformulation enabled us to prove the existence of stationary 
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R 




x 



Fig. 4. Spatial profiles of the amplitude R for fi = l,a = 1.5,(3 = — 1.2, L = 8.958 
from numerical simulation (dots) and the approximate solution (43) (solid line). 
The agreement is good, with the discrepancy mainly due to the long relaxation 
time close to the bifurcation. 

MAW solutions in the two limit cases corresponding to the bifurcation of the 
trivial solution A = (case I), or to the bifurcation of the plane wave solution 
of zero wavenumber (case II), when those solutions are within the Benjamin- 
Feir- Newell regime, or more generally in a region of instability denned by (21). 
That region coincides with the Eckhaus domain if \a\ = \/3\, but it is different 
otherwise. We proved the stability of MAW solutions for the full CGLe in a 
finite box with periodic boundary conditions in case II, where a homogeneous 
plane wave becomes unstable. We tested our analytical results by comparison 
of numerical integrations of the full CGLe with our approximate analytical 
solutions. 

In case I, unstable periodic hole solutions were shown to exist. This could 
not be inferred from any phase equation: around the defect point A = 0, the 
amplitude behaves non-analytically, namely piecewise affmely, and the phase 
is not defined. In case II we found the symmetric stable solutions observed 
in the numerical integrations of the CGLe just beyond the bifurcation point, 
using the box size L as the bifurcation parameter. This bifurcation was shown 
to be always supercritical in the Benjamin- Feir- Newell unstable regime. The 
MAWs continue to exist when the size L is increased. 

The analysis of MAWs bifurcating from a plane wave with wavenumber < 
q < 1 should be similar to the study of case II. It would be interesting to study 
the higher order instabilities of MAWs when the system size is increased be- 
yond the region in which our analysis takes place. It has been observed that 
stationary symmetrical MAWs bifurcate into uniformly-propagating asymmet- 
rical ones via a drift-pitchfork bifurcation. This happens when L is increased 
as a consequence of the growth of the amplitude of the modulation, and the 
increase of the spectral richness of the MAW solution. Moreover, MAWs are 
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expected to be the building blocks of phase turbulence, and the analytical 
analysis of their global stability may lead to a characterization of the sus- 
pected transition between phase and defect chaos in the CGLe [12,13]. 
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A Derivation of the governing equation 



We use (15) and (10) to rewrite (12) using only P and its spatial derivatives: 

P[aP xxx + bP x + cPP x ) x - P x (aP xxx + bP x + cPP x ) = (dP + e)P 3 (A.l) 

Note that this equation contains even numbers of derivatives of P in each term 
in parenthesis, and also that the powers of P increase while the derivatives 
decrease. We now rewrite the equation in a form which take advantage of this 
structure. For example, the following equation is equivalent to (A.l) for any 
real A: 

P(aP xxx + bP x + (c + \)PP X ) X - P x (aP xxx + bP x + (c + \)PP X ) 

= P\\P xx + dP 2 + eP), 

or, put in another form and introducing another real parameter k: 



(aP xx + bP + ^rP 2 ) x \ 9 a 

- I = XP XX + dP ■ eP ■ -kP , 



where we have written e + jk = e. 

In this equation, we have three free parameters: besides iv, introduced by the 
ansatz (4) as the carrier frequency of the solution, we have introduced free 
parameters A and k. We now fix A by imposing the condition 

a _ b _ c+A 

A"e" 2d ' { } 

which allows us to write the equation in a more suggestive form: 

the equation (17) that leads to the A-D ODE of section 3. 
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A is determined by (A. 2): 

A 2 + cA - 2ad = . (A.3) 

The discriminant of (A.3) is 
A = c 2 + 8ad 




1 + a 2 \ 2 ( 9(1 + a/3) 2 \ 
2((3-a)) [ ((3-a) 2 + ) ' 



So A > for any real values of a and (3, and the quadratic equation (A.3) 
always has two real roots 



. 3(1 + af3){\ + a 2 ) , 1 + a 2 r— — — — 

= 4(g-a)» 11 4GS^^ 9(1 + ^ +8( ^ -a) (A ' 4 ' 



Note that A is a function of a and (3 only. In some applications [19], the two 
values of A correspond to two distinct solutions of the CGLe. In our case, 
A is an intermediate variable used in the derivation and the proofs, but our 
solutions to the CGLe do not distinguish the two values of A. 



B Theorems used in the proofs 



We use successive approximation method to prove the existence of modulated 
amplitude waves. Below are listed several theorems from the theory of nonlin- 
ear oscillations taken from Hale's monograph [8]. 

Consider the system of equations 

z = Az + eZ{r,z,e) (B.l) 

where A is a constant matrix, e, r e R, and z, Z e M n . Z is a continuous 
function of r, z, e, periodic in r of period T. In the following, we only consider 
the case that Z is a smooth function. Without loss of generality, A can always 
be assumed to have the standard form 

A = diag(0 P)J B), 

Where P is a p x p zero matrix and B is a constant matrix with the property 
that the equation y = By has no nontrivial periodic solution of period T. 
Under these settings, if the successive approximation is applied to (B.l), we 
have 
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Theorem B.l Given d > b > 0, there is an e\ > such that for any given 
constant p vector a,\\a\\ < b and real e,\e\ < e i; there is a unique function 

z*{t) = z(t, a, e), with sup||2;*(r) || < d 

T 

which has continuous first derivative with respect to r and satisfies 

z * = Az* + eZ(r, z*, e) - eP Z(r, z*, e). 

Furthermore, z(r, a, 0) = a*, a* = col(a, 0), Po(z*) = a*, and z(r,a,e) has 
continuous first derivatives with respect to a, e. 

Po is defined as a projection operator on the Banach space S of continuous 
periodic functions of period T. If / e S, write / = col(g, h) where g is a p 
vector and h is a n — p vector, then 

P„(/)=col (t- 1 £g(t) dt,0j 

So, Po brings an element / in S to a constant vector which has the average 
values of g over one period as the first p components and zeros as the rest 
components. The equation satisfied by z* is different from (B.l) by a constant 
vector. By a proper choice of the starting vector a, we may make this con- 
stant vector zero to obtain a solution for the system (B.l). The mathematical 
statement is give by the following theorem. 

Theorem B.2 Let z(r, a, e) be the function given by the Theorem B.l for all 
\\a\\ < b < d, |e| < e\. If there exist an e 2 < e\ and a continuous function a(e) 
such that 

P Z(T,z{r,a{e),e),e) = 0, with \\a(e)\\ < b for \e\ < e 2 (B.2) 

then z(r,a(e),e) is a periodic solution of system (B.l) for ||e|| < e 2 . Con- 
versely, if system (B.l) has aperiodic solution z(r,e), of period T, \\z(t, e)|| < 
d, |e| < e 2; then z(t, e) = z(t, a(e), e). 

Therefore, the existence of a continuous function a(e) satisfying (B.2) is a 
necessary and sufficient condition for the existence of a periodic solution of 
system (B.l) of period T. As we do not know the exact functional form of the 
periodic solution, the condition (B.2) could not be solved explicitly. But by 
using implicit function theorem, we can show that the substitution into (B.2) 
of a proper approximate function of z(t, a, e) leads to the existence condition 
for periodic solutions. 

Theorem B.3 In the system (B.l), let 

Z = col{X,Y), z = col(x,y) 
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where X, x are p vectors and define 



1 r 1 

Xo(x, y,e) = - / X(t, x, y, e)d 



Jo 



It. 



If there is a p vector do, ||ao|| < d, such that 



X (a ,0,0) = 0, det 



dX (a ,0,Q) 



7^0 



(B.3) 



dx 



then there exists an e± > and a periodic function z{r,e), |e| < e\, of system 
(B.l) of period T with z(t,0) = col(a ,0). 

If we need to determine other parameters as a function of e in practical appli- 
cations, similar theorems could be derived. Specifically, in the main text we 
consider the period T as a function of e. It is clear that theorem B.3 applies 
if we suppose T(e) is continuous in e and bounded for |e| < e\. Furthermore, 
despite the use of the zeroth approximation in the above theorem, the nth 
approximation could be used instead. If simple (non-vanishing determinant) 
solutions to the determining equations can be found for e in the neighborhood 
of then system (B.l) has a periodic solution. 

If the system which we are studying possesses certain symmetries, we can 
prove the existence of particular symmetric solutions by a simplified version 
of determining equations. Let us define the symmetry first. 

Definition B.l Let z = f(r,z), where z,f'£ ~R n , be a system of differential 
equations. It is said to have the property E with respect to Q if there exists a 
nonsingular matrix Q such that 



where Pq is the projection operator defined before. 

Under this symmetry assumption the following theorems apply: 

Theorem B.4 Suppose Q = diag(Qi, Q 2 ) where Qi is a p x p matrix. If 
system (B.l) has property E with respect to this Q for all e. If a, \\a\\ < b, is a 
p vector and a* = col(a, 0) is a n vector, chosen in such a way that Qa* = a* , 
then the solution z(t, a, e) satisfies the relation 



Q 2 = I Qf(-r,Qz) = -f(r,z) QP = P Q 



Qz(-t, a, e) 



z(t, a, e) 



and consequently, 



Z(-r,z(-r,a,e),e) 



Qz(r,z(r,a,e),e) 
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Theorem B.5 If the j-th element of the diagonal of the matrix Q\ in Theo- 
rem B.4 is +1, then the j-th equation in the determining equations is equal to 
zero for every vector a* in Theorem B-4- 

The system (18) derived here from the 1-D CGLe has this symmetry, so the 
number of determining equations can be reduced using these two theorems. 
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